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We model the quasiparticle interference (QPI) pattern in the recently discovered (K,Tl)Fe2:Se2 superconduc- 
tors. We show in the superconducting state that, due to the absence of hole pockets at the Brillouin zone center, 
the quasiparticle scattering occurs around the momentum transfer q — (0, 0) and (±7r, ±7r) between elec- 
tron pockets located at the zone boundary. More importantly, although both d^2„y2-wave and s-wave pairing 
symmetry lead to nodeless quasiparticle excitations, distinct QPI features are predicted between both types of 
pairing symmetry. In the presence of a nonmagnetic impurity scattering, the QPI exhibits strongest scattering 
with q = (±7r, ±7r) for the ci^2_y2-wave pairing symmetry; while the strongest scattering exhibits a ring-like 
structure centered around both q = (0, 0) and (±7r, ±7r) for the isotropic s-wave pairing symmetry. A unique 
QPI pattern has also been predicted due to a local pair-potential-type impurity scattering. The significant con- 
trast in the QPI pattern between the cij.2_y2-wave and the isotropic s-wave pairing symmetry can be used to 
probe the pairing symmetry within the Fourier-transform STM technique. 
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Introduction. The very recent discovery of high-Tc (above 
30 K) superconductivity in AFe^Ses (A= K, Tl, Cs) llHl has 
generated new excitement in the condensed matter commu- 
nity. Relative to other iron-based superconductors (such as 
LaOFeAs, BaFe2As2, FeSe etc.), the end members, TlFe2Se2 
and KFe2Se2 (called as the "122" iron-selenides), are heav- 
ily electron doped (0.5 electron/Fe). Band structure calcu- 
lations iQ-H] on these end compounds show only electron 
pockets, primarily located around the M point of the Bril- 
louin zone (BZ) defined for a simple tetragonal structure. 
Angle-resolved photoemission spectroscopy (ARPES) mea- 
surements observed these electron-like Fermi surface pockets 
around the M points, and showed no hole-like pockets 
but very weak electron-like pockets ll iMlSll near the zone 
center F. The strong Fe-deficient compound (x < 1.6) shows 
insulating behavior IHIS. This is in contrast to other iron- 
based parent compounds, which are poor metals, raising the 
interest in the possibility of a Mott insulating state llSi - llSll in- 
duced by patterned Fe- vacancies \ 23). These observations 
add to the possibility that the pairing symmetry in the new 
compounds is unconventional. In particular, the absence of F- 
centered hole pockets would invalidate the popular s±-type of 
pairing symmetry, which was proposed for earlier iron-based 
superconductors. Recent calculations have predicted that the 
superconducting state could have dj,2_j,2-wave i2ll - l24ll and 
s-wave symmetry llSll . All these scenarios lead to nodeless 
superconducting gap structure, which is in agreement with the 
ARPES observations and other experiments. Because of the 
particular Fermi surfaces, conventional phase-sensitive mea- 
surements cannot be readily applied to differentiate the pair- 
ing states. 

Recently, one of the present authors and co-workers have 
proposed use of the existence or absence of intra-gap reso- 
nance states induced by a nonmagnetic impurity to probe the 
superconducting pairing symmetry 12511 . It has been found 
that the impurity-induced resonance state can only exist for 
a (i^2_j,2-wave pairing state. As mentioned above, since this 
dx^^y2-wa\e pairing state does not introduce nodal quasipar- 



ticles, due to the unique Fermi surface topology, the impurity- 
induced resonance state is located near the superconducting 
gap edge and requires a strong potential scattering. An alter- 
native technique, which can directly identify the sign change 
of the superconducting order parameter across various regions 
of Fermi surface, is the quasiparticle interference (QPI) probe. 
This technique has made a great stride in understanding low- 
energy quasiparticle properties and superconducting gap sym- 
metry in high-Tc cuprates ll26l - l28ll . The underlying principle 
for the QPI is that even a weak impurity scattering will mix 
two electronic states with two different momenta but on the 
same shell energy contour in the Brillouin zone, and the re- 
sultant momentum transfer (or scattering interference pattern) 
coiTesponds to the modulation in the local density of states, 
which can be measured by the Fourier-transform STM tech- 
nique ll29ll . The analysis of QPI in the presence of impurity 
scattering has been theoretically proposed Isol ? , 31 ] to probe 



the pairing symmetry in earlier iron-based superconductors. 
Results of later QPI measurements ||32ll are consistent with 
the scenario of the order parameter having a sign reversal 
across the electron and hole pockets. Recently, the QPI sig- 
natures have also been discussed either for the whole phase 
diagram including the metallic spin-density wave order llssll 
or in the presence of magnetic field 11341] . 

Motivated by this recent success, here we perform a de- 
tailed analysis of QPI in the newly discovered "122" iron- 
selenide superconductors. Both the nonmagnetic impurity 
scattering and pair potential scattering are considered. The 
latter type of scattering is more relevant to the STM exper- 
iments on samples with an applied magnetic field, in which 
Abrikosov vortices are generated. We show a pronounced 
difference in the QPI characteristic of a simple s-wave and 
dx2_y2-wa\e pairing symmetry. Because the Fermi surface of 
the new (K,Tl)Fe:cSe2 compounds comprise small pockets of 
only one type of carrier, this kind of study will also provide an 
opportunity to identify the unique role of Fermi surface topol- 
ogy in the QPI pattern even for the same pairing symmetry. 

QPI methodology. In view of the fact that there is either 
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no hole pocket or faint features of electron pockets in the su- 
perconducting "122" iron selenides, we consider here a sim- 
plified single-band model, which enables us to obtain a full 
understanding of the QPI pattern due to the different Fermi 
surface topology and pairing symmetry. Recently, a simi- 
lar simplified single-band model (Hubbard model) was also 
used to understand the magnetism in the iron-deficient com- 
pounds 1351]. The model Hamiltonian is defined on a two- 
dimensional (2D) square lattice and consists of the pristine 



and impurity scattering parts, H = Hq 
part Hq is. 
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Here the operators c^^ (ckcr) create (annihilate) an electron 
with momentum k and spin a. The quantity denotes the 
energy dispersion. We consider only the spin-singlet pairing 
and the superconducting gap function is described by Ak. The 
impurity scattering part is given by 

<T S 

(2) 

where the first term represents the zero-range normal potential 
scattering from a nonmagnetic impurity located at the origin 
of the lattice (for mathematical convenience) and the scatter- 
ing strength is denoted as Uq, while the second term describes 
the electron scattering due to a perturbation in the pair poten- 
tial of amplitude SA. For s-wave pairing symmetry, 775 = 1 
with the variable S in the summation taking only the value 
of zero; while for d-wave pairing symmetry ijs = 1 (^1) 
for 5 = zbi (±y). For our purpose here, we calculate the 
Green's function in the presence of either a nongmagnetic im- 
purity or pair potential impurity scattering, which is defined 
as Gfe-;T) = -(T,[^',(t)*]]) with *| = {4vC,l) as a 
two-component operator in the Nambu space. For the weak 
scattering potential, we make a first-order approximation and 
obtain the Green's function 

G{i,j;iuj„) = Go{i,j;ioJn) + ?7oG'o(i, 0; iw„)'r3Go(0, j; iaj„) 

+dA^risGo{i,Q;iu!n)TiGo{5,j;iuJri) , (3) 
5 

where a;„ = (2n + l)TrT is the Matsubara frequency with 
n an integer and T the electronic temperature, ri and T3 are 
the components of Pauli matrices in the Nambu space, while 
Go is the Green's function for the pristine system. From now 
on, we will term the usual nonmagnetic impurity scattering 
T3-scattering and the scattering off a pair potential impurity 
Ti -scattering. We emphasize again that the latter type of scat- 
tering is more relevant to the STM measurements of samples 
in the presence of vortices. The QPI is characterized by the 
Fourier-transform of the local density of states (LDOS), that 
is, pq{E) = Here the LDOS is given by 

Pi{uj) = — (2/7r)Im[(5Gii(i, i; iun E + i-y)], with Gn be- 
ing the site-diagonal normal (electron) component of the ma- 
trix Green's function. Note that we have measured the LDOS 



in the presence of the impurity scattering by removing the uni- 
form background. A little algebra yields 13611 : 

-[ME)Kk+n{E) + Kk{E)A+q{E)]} 

+ W ^ M[ME)Ki^+ci{E) + K^{E)A^+^{E)] 



-[B^iE)Ju+q{E) + ME)B^+^{E)]} 



(4) 



where the form factor i^k = 2 for s-wave pairing sym- 
metry, while Fk = 2(cosfcx — cos ky) for c?-wave pair- 
ing symmetry. The functions A, B, J, and K are de- 
fined as A{J) = -(2/7r)Im[Go.ii(i2)(k;F)] and B(A') = 
Re[Go^ii(i2)(k; E)]. Equation Q shows that the QPI pattern 
is determined by the convolution of the bare Green's func- 
tions in momentum space. For T3 -scattering, the convolution 
occurs between either normal Green's functions or anomalous 
ones (with a negative sign); while for ri -scattering, the con- 
volution is between normal and anomalous Green's functions. 
These bare Green's functions are given by 
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where the quasiparticle energy Fk — v ''k ' "k' 
v/(l + ek/i^k)/2 and = sign(Ak) - W^^k)/2 are 
the electron and hole parts of the Bogoliubov wavefunction 
amplitude. Therefore, for a given momentum k, the contri- 
bution to the Fourier amplitude Pq{E) is sensitive to the sign 
of Wk^k+q, which holds the key to reveal the uniqueness of a 
superconducting pairing symmetry in the QPI measurements. 

QPI in dj.2_y2 -wave and s-wave pairing symmetry. Before 
we present the numerical results on the QPI, we point out that, 
for unconventional pairing symmetry, the quasiparticle energy 
is very sensitive to the Fermi surface topology. To illustrate 
this point, we consider the normal state energy dispersion: 



^k = — 2t(cos kx + cos ky) — At' cos kx cos ky 



(6) 



where t and t' are the nearest-neighbor and next-nearest- 
neighbor hopping integrals and ji is the chemical potential. If, 
as relevant to high-Tc cuprates, we take t = l,t' = —0.3, and 
/i = —1.0, the Fermi pockets are centered around the (tt, tt) 
and equivalents in the BZ, as shown in Fig.fllal). These pock- 
ets are cut by the zone diagonals, which makes the quasipar- 
ticle excitations gapless in the dx2-y2-wa\e pairing symme- 
try, Ak = (Ao/2)(cosA:a; — cosky)/2, and a V-shaped pro- 
file of quasiparticle density of states around the Fermi energy 
(see Fig.[TJa2)). It enables the analysis of QPI in cuprates to 
provide a detailed band structure and dj.2_y2-wave gap struc- 
ture lieUii. However, if we take t = 0, t' = -1, and 
/i = —3, as a simplified modeling of the band structure in the 
"122" iron selenides, the Fermi pockets are centered around 
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FIG. 1. (Color online) Fermi surface contour and density of states 
witii t = 1, f' = -0.3, and = -1 (a) and t = Q, t' = -1, 
and /i = —3 (b), whicfi are relevant to the cuprates and "122" iron 
selenides, respectively. In panels (a2) and (b2), the solid lines are 
for d^2_y2-wave pairing symmetry (Ak = Ao(cos — cos ky)/2) 
with A() = 0.2), while the dashed lines are for isotropic s-wave 
pairing symmetry (Ak = Ao with Ao = 0.1). The dashed lines in 
panels (al) and (bl) represent the diagonals across the BZ. 

the (tt, 0) point and equivalents in the BZ (see Fig. [Itbl)). 
Due to the low density of electrons, these Fermi pockets are 
small in cross-section and, as revealed experimentally, are 
nearly isotropic. In this situation, the quasiparticle excitations 
are fully gapped and a well-shaped profile of density of states 
is obtained (see Fig. [T]i(b2)). For the isotropic s-wave pair- 
ing symmetry, the quasiparticle excitations are always fully 
gapped, irrespective of the detailed Fermi surface topology 
(see Fig. [na2)-(b2)). This explains why both d^2_y2-wave 
and isotropic s-wave pairing symmetry scenarios are compet- 
ing candidates for the photoemission spectroscopy measure- 
ments in "122" iron selenides. To clarify the pairing symme- 
try, the QPI would be a powerful technique. 

We now turn to the QPI by focusing on the special fea- 
tures arising from the dj.2_j,2-wave and isotropic s-wave pair- 
ing symmetries. We fix the band structure parameters with 
values of t = 0, t' = —1, and /i = —3, with the maxi- 
mum pair potential amplitude Aq ~ 0.1. A lattice size of 
Nl = 2048 X 2048 is typically taken. The strength of the 
T3- and Ti-scattering potentials are taken as Uq ~ 0.1 and 
SA = 0.01, respectively. Because for both d^2_y2-wave and 
isotropic s-wave pairing symmetry with the type of Fermi sur- 
face for "122" iron selenides, the quasiparticle excitations are 
fully gapped, we take the energy E ~ ±0.1, which is close 
to the coherent gap edge, for the QPI analysis. The intrinsic 
broadening parameter is taken as 7 = 0.02. 

In Fig. I2] we show the quasiparticle scattering pattern 
\pq{E) I in the presence of either rs-scattering or ri-scattering 
for the nodeless (1^2 _y2 -wave pairing symmetry. As can be 
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FIG. 2. (Color online) The quasiparticle scattering map, |pq(_E)|, 
for rg impurity scattering (upper panels) and ri impurity scattering 
(lower panels) at energy E — —0.1 (left column) and E — 0.1 (right 
column) with the nodeless d^2_y2-wave pairing symmetry. For the 
ra impurity scattering, the intensity is amplified by a factor of 10, 
while for the ri impurity scattering, it is amplified by a factor of 25. 

seen from Fig.|2](al)-(a2), for the ra-scattering, the QPI pat- 
tern exhibits a ring-like structure around momentum-transfer 
q = (0, 0) and q = (±7r, ±7r). However, the minimal inten- 
sity occurs at q = (0, 0), while the maximum intensity occurs 
at q = (±7r, ±7r). On the one hand, the maximum intensity 
at q = (±7r, ±7r) indicates strongest scattering between the 
Fermi pockets located at (±7r, 0) and (0, ±7r), which is due 
to the opposite signs of the superconducting order parameter 
around these two Fermi pockets. Furthermore, the ring-like 
structure centered around q = (±7r, ±7r) for —\E\ is much 
weaker in intensity than that for \E\. On the other hand, the 
ring-like structure around q = (0, 0) is contributed from the 
intra-pocket scattering. In the presence of ri -scattering, the 
quasiparticle scattering changes significantly when compared 
with the case of T3-scattering. In this case, the bright ring- 
structure only occurs around q = (0,0), while almost no mea- 
surable intensity is obtained around q = (±7r, ±7r). In partic- 
ular, the maximum intensity is located at the point q = (0,0). 
The significant difference in the quasiparticle scattering pat- 
tern between the cases of rs-scattering and ri-scattering sug- 
gests that the QPI is also very sensitive to the nature of the 
impurity scattering. Therefore, in the interpretation of exper- 
imental data for superconducting pairing symmetry, caution 
must be taken as to whether the impurity scattering is of T3 or 
Ti nature. 

In Fig. [3] we show the quasiparticle scattering pattern 
Pq(i?) I in the presence of either r3 -scattering or ti -scattering 
for isotropic s-wave pairing symmetry. For T3 -scattering (see 
Fig. [3lal)-(a2)), the ring-like structure appears around both 
q = (0, 0) and (±7r, ±7r). However, in striking contrast to 
the case of dx2^y2-wave pairing symmetry (in the same con- 
dition of Ts -scattering), the maximum intensity is located di- 
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FIG. 3. (Color online) The QPI maps for the rg impurity scattering 
(upper panels) and ri impurity scattering (lower panels) at energy 
E = —0.1 (left column) and E = 0.1 (right column) with isotropic 
s-wave pairing symmetry. For rs impurity scattering, the intensity 
is amplified by a factor of 10, while for ri impurity scattering it is 
amplified by a factor of 50. 



pairing symmetry due to the sign change in the superconduct- 
ing gap across electron pockets at the BZ boundary; while 
the strongest scattering exhibits a ring-like structure centered 
around both q = (0, 0) and q = (±7r, ±7r) for isotropic s- 
wave pairing symmetry. In the presence of ti impurity scat- 
tering, the strongest QPI intensity occurs only at q = (0, 0) 
for the dx2_y2-wa\e pairing symmetry, while it occurs at both 
q = (0, 0) and (±7r, ±7r) for isotropic s-wave pairing sym- 
metry. This analysis shows the sensitivity of the QPI pattern 
to the nature of impurity scattering. The significant contrast 
in the QPI pattern between the d^2_y2-wave and the isotropic 
s-wave pairing symmetry in the presence of the same type 
of impurity scattering can be very efficient for probing the 
pairing symmetry in the "122" iron selenide superconductors 
within the Fourier-transform STM technique. 
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LDRD-DR Program, and the DOE Office of Basic Energy of 
Sciences. 



rectly on the rings, while almost no visible intensity is ob- 
tained at both q = (0, 0) and (±7r, ±7r) points. For ti- 
scattering (see Fig.[3jbl)-(b2)), the ring-like structure remains 
centered around the q = (0,0) and q = (±7r, ±7r) points. 
More noticeably, the maximum intensity is now located at 
both q = (0,0) and (±7r, ±7r) points, which is also signifi- 
cantly different than the case of d^2_y2-wave pairing symme- 
try, where the maximum intensity occurs only at q = (0,0). 
Also the difference in the QPI pattern between rs -scattering 
and Ti -scattering indicates the importance of identifying the 
nature of impurity scattering when QPI data are interpreted to 
probe the superconducting pairing symmetry. 

Finally we note that the QPI analysis for the earlier 
iron pnictide superconductors with s±-wave pairing symme- 
try OOl [31I [33fl shows the most pronounced scattering at 
q = (±7r, 0) [(0, ±7r)], which is between the electron (at the 
M point in the BZ) and hole (at the F point in the BZ) pock- 
ets. For the "122" iron selenide superconductors, due to the 
absence of the hole pockets at the BZ center, the quasiparti- 
cle scattering can only occur between the electron pockets and 
the QPI pattern shows pronounced structure around q = (0,0) 
[(±7r, ±7r)] for either d^2_j,2-wave or isotropic s-wave pair- 
ing symmetry. 

Conclusion. We have studied the quasiparticle interference 
pattern due to both T3 and ri impurity scattering in the re- 
cently discovered "122" iron selenide superconductors. We 
have shown in the superconducting state that, although both 
dj.2_y2-wa\e and s-wave pairing symmetry lead to nodeless 
quasiparticle excitations, the QPI pattern is strikingly different 
between the two types of pairing symmetry. In the presence 
of Ts impurity scattering, the QPI exhibits strongest scatter- 
ing with momentum transfer q = (±7r, ±7r) for rfj.2_,y2-wave 
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